The effect of an in-plane magnetic field on the interlayer transport of quasiparticles in 

layered superconductors 
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We consider the quasiparticle c-axis conductivity in highly anisotropic layered compounds in 
the presence of the magnetic field parallel to the layers. We show that at low temperatures the 
quasiparticle interlayer conductivity depends strongly on the orientation of the in-plane magnetic 
field if the excitation gap has nodes on the Fermi surface. Thus measurements of the angle-dependent 
c-axis (out-of-plane) magnetoresistance, as a function of the orientation of the magnetic field in the 
layers, provide information on the momentum dependence of the superconducting gap (or pseudogap) 
on the Fermi surface. Clean and highly anisotropic layered superconductors seem to be the best 
candidates for probing the existence and location of the nodes on the Fermi surface. 



PACS: 74.25.Fy, 74.70.Kn, 74.72.-h 

The symmetry of the excitation gap in the supercon- 
ducting and normal state of the cuprates, the quasi-two- 
dimensional organic salts, and the ruthenate oxides has 
been the focus of theoretical and experimental studies 
for the last several years. Now there is consensus that 
most of the cuprates are not only anomalous metals that 
develop a pseudogap when underdoped m , but also that 
they are unconventional (d-wave) superconductors with 
nodes on the Fermi surface, as demonstrated in phase- 
sensitive Josephson junction experiments J2|. The situa- 
tion is less clear in the organic salts and the ruthenates, 
in spite of many reports of power laws in the tempera- 
ture behavior of various transport and thermodynamic 
properties (see, e.g., Refs. |],§j]). Spin fluctuation models 
for the organics predict a d xy superconducting state simi- 
lar to the cuprates ||||] . Since phase-sensitive Josephson 
junction and angle-resolved photoemission spectroscopy 
experiments are not available for these materials, other 
more stringent experiments are required. 

In this Letter we propose angle-dependent magnetore- 
sistance oscillation (AMRO) experiments which directly 
probe the locations of the line nodes of the gap in the 
quasiparticle excitation spectrum in the superconducting 
phase or in the normal state with a pseudogap. We show 
that in highly anisotropic layered materials the depen- 
dence of the quasiparticle c-axis PV characteristics on the 
orientation of a magnetic field parallel to the layers, Bii, 
allows one to extract the momentum dependence of the 
gap function. In particular, we calculate the quasiparti- 
cle c-axis conductivity, a q , in Josephson coupled layered 
gapless superconductors (e.g., rf-wave or p-wave pairing) 
as a function of By = (B x , B y ) and compare it with re- 
sults for s-wave pairing. We assume a significant contri- 
bution of coherent electron tunneling between adjacent 
layers. The condition for coherent tunneling in the pres- 
ence of an in-plane magnetic field requires |<5k| <C |b|, 
where |<Sk| is the change of the in-plane momentum of a 
tunneling electron, where b = (2ns/^o)(B y ,—B x ), and 



s is the interlayer spacing. We show that under such 
conditions the dependence of a q on the orientation of b 
enables one to determine whether the superconducting 
gap (or any other kind of gap) has nodes on the Fermi 
surface, as well as information on the location of these 
nodes. The conditions that are necessary to apply our 
method seem to be fulfilled in highly anisotropic and 
very clean organic salts like the quasi-two-dimensional 
BEDT-TTF [bis(ethylenedithio)-tetrathiafulvalene] su- 
perconductors B , the Bi- or Tl-based cuprates MM , and 
the oxide superconductor Sr 2 Ru04 ||]. 

For perfect crystals with translational invarian.ee we 
use the interlayer Hamiltonian, which describes the con- 
servation of the in-plane momentum when electrons tun- 
nel between the layers. In Josephson coupled supercon- 
ducting layers the in-plane magnetic field penetrates al- 
most freely into the sample, inducing the vector potential 
a z (r) = b • r, where r = (a;, y) is the in-plane coordinate. 
Thus the interlayer Hamiltonian or tunneling Hamilto- 
nian can be written as 



H± = t A 
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where ip n a is the annihilation operator for electrons in 
layer n with spin a, and t ± is the interlayer transfer in- 
tegral. We assume that t ± is isotropic in the layers. The 
vector potential a z (r) leads to a change in the in-plane 
momentum of the tunneling electron, k — > k + b, where 
k = (k x ,ky) is the in-plane momentum. In the momen- 
tum representation the interlayer tunneling Hamiltonian 
has the form 



H±(b) = t ± ^2 KU<xVVK,k+b,<x + h.c] . 



(2) 



Defects in the crystal structure lead to spatial variations 
in t±. The spatial average, (£j_), determines the coher- 



ent tunneling part. We assume that it dominates the 
low-temperature c-axis transport and use in (1) and (2) 
t± = (t±), for details see Ref. g. Scattering inside the 
layers leads to a change of the momentum while tun- 
neling, <5k| = 1/^||. Here In is the effective mean-free- 
path for scattering of quasiparticles inside the layers. In 
conductivity measurements a finite voltage V is applied 
between neighboring layers and hence the energy of the 
tunneling electron changes by V. In the limit V — > the 
energy is conserved within the accuracy of the tempera- 
ture T . The important point is that at low temperatures 
only quasiparticles near the nodes of the gap on the Fermi 
surface can contribute to the dissipative quasiparticle c- 
axis transport. Let us denote the positions of the nodes 
on the Fermi surface by the momenta k g . Then at low 
T and b 3> l/£y the contribution of the quasiparticles 
near the nodes to the conductivity is significant when 
the change of the quasiparticle energy, -Ek+b — E^, is 
small, i.e., when hvf(k g ) ■ b <C T <C A , see Fig. |[ Here 
k = k/fc is a unit vector, v/ is the Fermi velocity, and Ao 
is the amplitude of the energy gap. This means that the 
dominant contribution to a q comes from quasiparticles 
in the vicinity of the nodes and contribution of a given 
node kg to a q is maximal for b _L v/(k s ). 

In the superconducting state the Josephson current 
contributes to the c-axis transport. Thus to obtain only 
the quasiparticle contribution the Cooper pair current 
must be suppressed when the interlayer current is mea- 
sured as a function of the orientation of the in-plane 
magnetic field. The suppression of the Josephson cur- 
rent can be achieved by applying a c-axis current signifi- 
cantly exceeding the Josephson critical current M . Such 
a condition is easily fulfilled, considering that the in-plane 
magnetic field strongly diminishes the Josephson critical 
current, without significantly affecting the superconduc- 
tivity inside the layers. Another method for observing 
the quasiparticle current is to perform the measurements 
of the c-axis conductivity in the resistive state of the in- 
trinsic Josephson junctions |10[ . 

We calculate the quasiparticle conductivity a q using 
the BCS theory for d-wave and p-wave pairing in the 
presence of elastic, isotropic in-plane scattering, char- 
acterized by the normal-state scattering rate r <C Ao. 
Here Ao is the amplitude of the superconducting spin- 
singlet d-wave gap [A(k) = A (T)(fc 2 - k%) or A(k) = 
2Ao(T)k x k v ] or of the spin-triplet p-wave gap along 
a given spin direction [A(k) = Ao(T)k x or A(k) = 
Ao(T)k y }. We consider first a circular Fermi surface in- 
side the layers and a cylindrical, three-dimensional open 
Fermi surface with t ± much smaller than any other rel- 
evant energy scale in the system. For simplicity, we ne- 
glect the effect of the Zeeman term on the in-plane Green 
functions uM. Although it will change the field depen- 
dence of a q , it will not affect its angular dependence. 
The perturbation theory with respect to H± gives the 



expression for the quasiparticle conductivity in terms of 
the quasiparticle spectral functions |12| : 
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where ImG(k, lu) = — 7rA(k, lu) is the imaginary part of 
the Green function. Within the BCS theory we obtain 



A(k,w) 



(i + &/£kh 



(i - a/^kh 
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where E k = [££ + A 2 (k)] 1 / 2 , and £k is the quasipar- 
ticle energy in the normal state, measured from the 
Fermi level. The scattering rate of the quasiparticles is 
7(c<j) = — 2ImS(o;), where £(w) is the diagonal part pro- 
portional to the unit matrix of the impurity self-energy, 
and lu'(lu) = lu — ReE(oj). The impurity self-energy has 
to be calculated self-consistently and is energy (temper- 
ature) dependent, even for elastic scattering, in the su- 
perconducting state |13| ]. 

In unconventional superconductors even nonmagnetic 
impurities lead to gapless states in the vicinity of the 
nodes of the gap function on the Fermi surface. At low 
temperatures, T <C 7(0) <C A o (0), the dominant contri- 
bution to the quasiparticle conductivity comes from these 
gapless states. In rather clean superconductors and at 
low energies 7(0;) becomes 7(0) ~ [rAo(O)] 1 / 2 in the case 
of strong scattering, and 7(0) ~ A o (0) exp[— A o (0)/r] 
for weak scattering (Born limit). Note that the relevant 
phase space for scattering of quasiparticles is restricted 
to a fraction ~ 7(0)/Ao(0) of phase space. At low tem- 
peratures, T <ti 7(0), we obtain 



0g(b)«^-^ fd 2 kA{k + b,0)A(k,0). 
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In the following, we take into account that regions 
near the nodes contribute mainly to the integral over 
k. We replace E 2 z+h w £ 2 +b + A 2 (k), where £ k +b ~ 
£k + hh ■ Vf (k) , and parameterize k = (cos cp, sin <p) and 
b = b(cos9, sin^). For example, the N g = 4 nodal angles 
for d K 2_ a 2 pairing are given by tp g = (2g— l)7r/4. Finally, 
we find at low temperatures and at low fields a universal 
conductivity, similar to the in-plane transport uM , 
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(6) 



with a g = (£n/2vf)h-Vf(is. g ) = (tnb/2)cos(cp g —6), where 
£» = h,Vf/'y(0) is the effective scattering mean-free-path, 
and Jo = 2et 2 L Nf/h is the Josephson critical current den- 
sity. Here Nf = m/(2ir 2 h 2 ) is the 2D density of states 



per spin of quasiparticles with mass m. This result is 
readily generalized to an elliptical 2D Fermi surface by 
rescaling the momenta fe, = ^Jm^jmk'^ (i — X,y), where 



are the effective masses and 



m^m y . Then 



Eq. (pi) is valid with a g — (£u 6/2) W/m/m x cos 6 cos (p g + 
yjvfijm y sin 9 sin ip g \ . 

The largest contribution of a given node to the con- 
ductivity comes when b _L vj(k 9 ). The conductivity is 
maximal for b along the nodes with periodicity 9 = 7r/2, 
as shown in Fig. g For low fields, (l\\b/2) 2 < 1), a q (h) 
shows a quadratic field dependence. In the high field 
limit, (£||6/2) 2 3> f , the conductivity falls off as b~ 2 In b. 

In the case of 7(0) <C T <C Ao(0), the quasiparticle 
scattering rate is energy (temperature) dependent and 
7(T) ~ r(T/A ) n with n = 1 in the weak scattering 
regime (Born limit), and n = — 1 for strong scattering 
||13|| . Thus for fields b <C Aq(0) /hv/ we obtain after inte- 
grating over to 
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In the next step, we expand -E k +b — £"k ~ (OE^/dk) ■ b = 
E^ 1 £khv f (k) ■ b, and take into account that at low tem- 
peratures, T <C Ao(0), only small values of £ k and angles 
near the nodes, (p g , contribute to the quasiparticle con- 
ductivity. Then we obtain 
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(8) 



where now ^|| ~ vjAq/TT in the Born limit and lu ~ 
v/T/rAo in the strong scattering limit. In the high field 
limit, (£||6/2) 2 ^> 1, the quasiparticle c-axis conductivity 
falls off as 1/6. 

These results are quite general and also valid for p- 
wave superconductors. In this case tp g = 0, 7r or ip g = 
7r/2, 37r/2. Now a q (6) is minimal for fields along the 
nodes with periodicity w. A p-wave state has larger oscil- 
lations than a (i-wave state for the same parameter £116. 

Angle-dependent oscillations of a q (9) were observed in 
the normal state in organic conductors [[L5| , cuprates M , 
and ruthenates [p[ , and were explained in terms of Fermi 
surface anisotropy. In materials with an elliptical 2D 
Fermi surface the angular periodicity of <J q [0) is n. The 
amplitude of the oscillations is determined by the param- 
eter Svfb/T, where Svf is the variation of the anisotropic 
Fermi velocity. 

The periodicity and height of the AMRO, induced 
by the topology of the Fermi surface in the normal 
state, changes when a superconducting gap or pseudogap 
emerges. If an orthorhombic (organic) crystal undergoes 
a superconducting transition with s-wave pairing, the 
7r-periodicity will preserve, but <j q will tend to zero as 
T — > 0. For a d-wave state an additional 7r/2-periodicity 



emerges on cooling and a q saturates as T — > 0. In tetrag- 
onal cuprate crystals the periodicity is tt/2, both in the 
superconducting and normal state M. The positions of 
the AMRO maxima in the normal and superconducting 
states may differ (as in Tl-2201 fiffl]) or coincide. In the 
former case, the maxima observed in the normal state 
should drop and new ones should appear on cooling be- 
low T c and further increase with increasing field b. In the 
latter case, the positions of the AMRO maxima do not 
change below T c , but a q saturates on cooling at T — > 
for d-wave superconductivity, while for s-wave pairing u q 
vanishes exponentially. The main parameter which deter- 
mines the angular behavior of cr q (8) at low temperatures, 
T < T c , is £\\b = Vfb/-f. We note that 7 > T, while 5v f 
can be of order Vf. Thus, we anticipate that the am- 
plitude of oscillations will be about the same above and 
below T c in tetragonal crystals and weaker in orthorhom- 
bic crystals. 

Using typical values for the interlayer spacing, s rs 
15 A, in organic and cuprate superconductors, e.g., 
Bi2Sr2CaCu208-5 (Bi-2212) crystals, we obtain b w 
4.56 B (T fxm) ~ l . Then at low temperatures a mean-free- 
path of £11 > 200 A is needed for a d-wave state in order to 
observe significant oscillations, i.e., oscillations that are 
larger than 1% of the total magnetoresistivity for mag- 
netic fields of order 10 T. It is feasible that this condition 
is fulfilled in the organic BEDT-TTF superconductors, 
where £11 ~ 1000 A was deduced from de Haas- van Alphen 
and Shubnikov-de Haas measurements Eq]. In the high- 
temperature superconductors Bi-2201 and Bi-2212 the 
mean-free-path, t\\ , of quasiparticles at low temperatures 
is sufficiently long Mm to observe this effect, too. 

Until now we have discussed the dependence of the 
quasiparticle conductivity on the orientation of b. More 
generally, the dependence of a q on b, described by 
Eq- @; provides information on the spatial correlations 
in the system. Namely, if cr ? (b) is obtained experimen- 
tally then the spatial dependence of the spectral density 
of the Green function may be extracted by an inverse 
Fourier transformation with respect to b, 

r + °° u 

(1/4T) / du cosh" 2 — A(r, u>)A(-r, u) = 



2T 
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In the low temperature regime the left-hand side is sim- 
ply A(r,0)A(-r,0) = A 2 (r, 0) and the coordinate de- 
pendence of the spectral function at zero energy may 
be calculated directly from measurements of cr q (b). A 
similar method was used previously to obtain informa- 
tion on the correlation function of pancake vortices in 
a vortex liquid from the dependence of the plasma res- 
onance frequency or from the dependence of a± on the 
magnetic field parallel to the layers with a fixed out-of- 
plane component p^Jlq]. For superconductors, in the 



absence of a perpendicular magnetic field, the coordi- 
nate dependence of A(r, 0) provides information on the 
superconducting correlation length. At low temperatures 
in highly anisotropic layered metals in the normal state 
the combination of parallel and perpendicular magnetic 
field components leads to oscillations of the c-axis con- 
ductivity as a function of & or as a function of the angle 
between the magnetic field and the c-axis Jl9|,^0| . 

The coordinate dependence of the spectral density of 
the Green function at arbitrary energy may be extracted 
from c-axis I-V measurements in the presence of an in- 
plane magnetic field. The c-axis current density is 



J g (b,V) 
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In the limit of low temperatures, T — > 0, the spectral den- 
sities are only slowly varying with respect to u>. In this 
limit we solve the integral equation (nQ) for A(r,u>) at 
known J q (h, V) by Fourier and Laplace transformation, 



/>oo 
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Jo 
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In conclusion, we have shown that, in principle, at low 
temperatures the spatial dependence of the spectral func- 
tion A(r, to) can be determined from measurements of 
the c-axis I-V characteristic as a function of the mag- 
netic field Bii. The dependence of a q on the orientation 
of the in-plane field Bm provides information on the ex- 
istence of line nodes on the Fermi surface as well as on 
their locations. Measurements of the angle-dependent 
in-plane magnetoresistance will be useful to map out the 
momentum distribution of the superconducting gap or 
pseudogap in organic salts and cuprates. 

We thank V.M. Yakovenko and A.V. Balatsky for 
many useful discussions. This work was supported by 
the Los Alamos National Laboratory under the auspices 
of the U.S. Department of Energy. 



[1] M.R. Norman et al, Nature 392, 157 (1998). 

[2] C.C. Tsuei and J.R. Kirtley, Physica C 282-287, 4 
(1997); C.C. Tsuei et al, Nature 387, 481 (1997). 

[3] K. Kanoda, Hyperfine Interact. 104, 235 (1997). 

[4] T. Ishiguro, K. Yamaji, and G. Saito, Organic Supercon- 
ductors, 2nd ed., Springer- Verlag, Berlin (1998). 

[5] J. Schmalian, Phys. Rev. Let t. 81, 4232 (1998). 

[6] Yu.I. Latyshev et al, e- print 3ond-mat/9903256. 

[7] N.E. Hussey et al, Phys. Rev. Lett. 76, 122 (1996). 

[8] E. Ohmichi et al, Phys. Rev. B 59, 7263 (1999). 

[9] M. Suzuki, S. Karimoto, and K. Namekawa, J. Phys. Soc. 
Jpn. 67, 732 (1998). 



[10] 

[11] 
[12] 

[13] 



[14] 



[15] 

[16] 
[17] 

[18] 
[19] 

[20] 



R. Kleiner and P. Miiller, Phys. Rev. B 49, 1327 (1994). 

K. Yang and S.L. Sondhi, Phys. Rev. B 57, 8566 (1998). 

G.D. Mahan, Many-Particle Physics, 2nd ed., chap. 9, 

Plenum Press, New York (1990). 

C.J. Pethick and D. Pines, Phys. Rev. Lett. 57, 118 

(1986); P.J. Hirschfeld, D. Vollhardt, and P. Woffle, Sol. 

State Comm. 59, 111 (1986). 

PA. Lee, Phys. Rev. Lett. 71, 1887 (1993); M.J. Graf, 

S.-K. Yip, J. A. Sauls, and D. Rainer, Phys. Rev. B 53, 

15147 (1996). 

A.G. Lebed and N.N. Bagmet, Phys. Rev. B 55, R8654 

(1997). 

N. Toyota et al, Sol. State. Comm. 72, 859 (1989). 

A.E. Koshelev, L.N. Bulaevskii, and M.P. Maley, Phys. 

Rev. Lett. 81, 902 (1998). 

N. Morozov et al, Phys. Rev. Lett. 82, 1008 (1999). 

R.H. McKenzi e and P. Moses, Phy s. Rev. Lett. 81, 4492 

(1998); e-print |cond-mat/9812il3[ 

A. Dragulcscu, V.M . Yakovenko, and D.J. Singh, e-print 



cond-mat/9811101 




FIG. 1. Sketch of the shifted Fermi circles in adjacent 
layers, which a tunneling electron sees, due to an applied 
in-plane magnetic field b. The filled circles indicate the 
nodal regions of a d xy order parameter with nodal vectors 
|ki|«|ki|and|k3| = |ki|-|b|. 
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FIG. 2. AMRO, p = 1/o-q, at low temperatures, 

T <C 7(0), for a d a .2_ H 2-wave superconductor at different pa- 
rameters l\\ b = 1, 2,4 and for an elliptical Fermi surface with 
specified mass anisotropies. 



